ALMOST COMMUTING MATRICES, LOCALIZED WANNIER 
FUNCTIONS, AND THE QUANTUM HALL EFFECT 



MATTHEW B. HASTINGS AND TERRY A. LORING 



Abstract. For models of non-interacting fermions moving within sites arranged on a sur- 
face in three dimensional space, there can be obstructions to finding localized Wannier 
functions. We show that such obstructions are if-theoretic obstructions to approximating 
almost commuting, complex-valued matrices by commuting matrices, and we demonstrate 
numerically the presence of this obstruction for a lattice model of the quantum Hall effect in 
a spherical geometry. The numerical calculation of the obstruction is straightforward, and 
does not require translational invariance or introducing a flux torus. 

We further show that there is a index obstruction to approximating almost com- 
muting self-dual matrices by exactly commuting self-dual matrices, and present additional 
conjectures regarding the approximation of almost commuting real and self-dual matrices by 
exactly commuting real and self-dual matrices. The motivation for considering this problem 
is the case of physical systems with additional antiunitary symmetries such as time reversal 
or particle-hole conjugation. 

Finally, in the case of the sphere — mathematically speaking three almost commuting 
Hermitians whose sum of square is near the identity — we give the first quantitative result 
showing this index is the only obstruction to finding commuting approximations. We review 
the known non-quantitative results for the torus. 
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Almost Commuting Matrices 



1. Asymptotic Commutants of Finite Rank Projections 

Given a list of bounded operators on infinite dimensional Hilbert space, it is often natural 
to seek a finite rank projection P that almost commutes with that set. The C*-algebraist 
would do so in the study of quasidiagonality, O [16], HQ]. In physics, we are interested in a 
projection onto a band of energy states separated from the rest of the spectrum by an energy 
gap; assuming the underlying Hamiltonian is local, this projection will itself be local due to 
the gap, and hence will approximately commute with a list of observables. 

Whatever exact relations might be known to hold for the original operators (Hi, . . . , H r ) 
will generally hold only approximately for the compressions (PX\P, . . . , PX r P). In a lattice 
model, the projection might be from a finite dimensional space to a space whose dimension 
is much lower, but still the outcome is finite-rank operators that approximately satisfy some 
relations. Can these be approximated by finite-rank operators that exactly satisfy those 
relations.? 

For example, if Xi and X 2 in B(H) satisfy —I < Xj < I and [Xi, X 2 ] = 0, then P almost 
commuting with the Xj implies 

-I < PXjP < I, 
\\[PX 1 P,PX 2 P)\\ «0. 

We especially want to know if these almost commuting Hermitian operators are close to 
commuting Hermitian operators in the corner PM(M)P = Mfc(C). It is sufficient to answer 
this question: can two almost commuting Hermitian matrices be approximated by commuting 
Hermitian matrices? The answer is yes. This is Lin's theorem [30] . 

The situation very different if we consider three almost commuting Hermitians. Specifi- 
cally, Theorem 13. 181 shows that there are X, Y, Z in B(H) satisfying 

[X,Y} = [X,Z] = [Y,Z}=0, 
X 2 + Y 2 + Z 2 = I 

and finite rank projections P n asymptotically commuting with X, Y, Z and yet such that 
there do not exist triples (H[, H 2 , H' 3 ) of operators with 

[H' r , H' s ] = 0, 
(H'if + (H' 2 f + (H' 3 f = I 
\\H r -H' r \\ «0, 

where Hi = PXP, H 2 = PYP, and H% = PZP. We will make this precise below, but the 
example is a variation on the examples in [8J (33J [39] . 

The key to showing this result is the presence of an index obstruction. Conversely, our main 
quantitative result is Theorem 13. 16[ which gives quantitative bounds on how accurately three 
almost commuting Hermitian matrices with vanishing index obstruction can be approximated 
by exactly commuting matrices. Specifically, we show that 

Theorem 1.1. Suppose (Hi, H 2 , H3) is a 5 -representation of the sphere by matrices (defined 
below). If the index bott(ifi, H 2 , H 3 ), defined below, is vanishing, then there are commuting 
Hermitian matrices H[,H 2 ,H' 3 with 

(H'if + (H' 2 f + (Htf = I 



Almost Commuting Matrices 



3 



and 

\\H' r -H r \\ < e{5) 

for all r, where e(8) = E(1/S)S 1 ^ 12 and the function E(x) grows more slowly than any power 
of x. 

Unless stated otherwise, matrices and vector spaces are over the complex numbers. We 
use here X T to denote transpose and X* the conjugate transpose. We always use || || to 
mean the operator norm. Operator means a bounded linear operator on C n or Hilbert space. 
Contraction means a operator of norm at most one. 

Example 1.2. We give an example of n- dimensional triples of matrices which almost com- 
mute but are not close to exactly commuting matrices. Let n — 25 + 1, where 5* is either 
integer or half-integer. 

Consider the spin matrices S^S^S* 3 for a quantum spin 5 with Hi = 5 1 / a/ 5(5 + 1), 
H 2 = S 2 /^S(S + 1), H 3 = S 3 /^/S(S + l). Thus, 



S/y/S(S+l) 



(5-1)^5(5 + 1) 



-S+l)/y/S(S+l) 



-S/y/S(S + 1) . 



and 

[H r ,H s ] = e rst H t / v /S(S + l), 

where e rst is a totally anti-symmetric tensor with e 123 = 1. For a mathematically oriented 
reader, the 5 matrices are a representation of the Lie algebra of 577(2). 
Note that 

Hi + Hi + Hi = I. 

It is easy to see that 

||i7 r , H s \\ < 1/5. 

Lemmas 13.51 and 13.81 combine to tell us that if H[, H' 2 and H' z are commuting, n-by-n 
Hermitian matrices then 



\\H[ - + \\H^ - H 2 \\ + \\H' 3 - H 3 \\ > VI "4/5. 

Choi [7j produced a slightly better estimate with essentially the same matrices. 

These spin matrices form an example of what we call an approximate representation of 
the sphere. The mathematics oriented reader should think about generators and relations 
for the C*-algebra Cq(S 2 ). The physics oriented reader should think that we describe the 
coordinates of a particle moving on the surface of a sphere; in the presence of a magnetic 
field, the particle's position is blurred out on the scale of a magnetic length and the different 
coordinates cease to commute exactly. 

To all manner of surfaces in M d , or C d , there are associated C*-algebras and related collec- 
tions of almost commuting d-tuples of matrices. We focus on the surfaces most prominent 
in physical models: the disk, square, annulus, cylinder, sphere and torus. 
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Definition 1.3. Suppose 5 > 0. A triple (Hi, H 2 , H 3 ) of operators or matrices is called a 
5 -representation of the sphere if 

H: = H r , (Vr) 
|| [H r ,H s ] || <5, (W^s) 
\\Hl + Hi + Hj - I\\ < 5. 

Definition 1.4. Suppose 5 > 0. A pair (Ui, U 2 ) of operators or matrices is called a 5- 

representation of the torus if 

U* r U r = U r U* = I, (Vr) 
\\[Ui,U 2 ]\\<6, 

Definition 1.5. Suppose 5 > 0. A pair (Hi,H 2 ) of operators or matrices is called a 5- 

representation of the square if 

-I <H r < I, (Vr) 
II [Hi,H 2 ] || < 5, 

Definition 1.6. Suppose 5 > 0. An operator or matrix X is called a 5 -representation of the 
disk if 

||X|| < 1, 
\\[X*,X]\\<5, 

Definition 1.7. Suppose 5 > 0. An operator or matrix X is called a 5 -representation of the 
annulus if 

-/ < \X\ < I, 
2 - I I - > 

||[X*,X]||<5, 

Definition 1.8. Suppose 5 > 0. A pair (U,K) of operators or matrices is called a 5- 

representation of the cylinder if 

[/*[/ = J/ET = /, 
-I <K <I, 
\\%K}\\<5, 

In the above, we will usually say "exact representation" instead of "0-representation." If 
we have a ^-representation for 5 > and don't wish to emphasize the exact value of 5, we 
will say "approximate representation." 

We will define an invariant, called the Bott index, that applies to approximate representa- 
tions of the sphere. This is an invariant that distinguishes those that can be approximated 
by exact representations and those that cannot. 

There is a more general index, defined explicitly using ii'-theory, that applies to almost 
commuting triples in C*-algebras. This has been studied in many papers, including |4^[TT|[3T]. 
Where possible we offer direct proofs in the language of matrix theory, with careful error 
estimates and avoiding .fT-theory or C*-algebras. 

The Bott index was discovered first in the context of ^-representations of the torus, i.e. 
for almost commuting unitaries U and V. We give four descriptions of this index. Their 
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complexity varies, but the simplest fits in one sentence. If VU is close to UV the determinant 
applied to a short path between UV and VU will create a closed path in the punctured plane 
whose winding number equals the Bott index of (U, V). 

The paper is organized as follows: in the next two sections we prove results needed for 
Theorem 13.161 In section 4 we review non-quantitative results on the torus. In section 5 
we consider physics applications, and in section 6 we consider index obstructions to approx- 
imation of almost commuting real and self-dual matrices by exactly commuting real and 
self-dual matrices, and we present a Z 2 obstruction in the self-dual case. 

2. Matrices that Almost Represent the Disk or Annulus 

There is no second cohomology for the disk, square, annulus or cylinder. This means there 
will be no obstruction (other than hard work) to perturbing approximate representations to 
exact representations. We are able to make this precise in quantitative theorems. 

We start with a minor variation to the quantitative version of Lin's theorem in [18J. 



Theorem 2.1. Suppose (H\,H 2 ) is a 5 -representation of the square by matrices. Then, 
there exists an exact representation (Kx,K 2 ) of the square with 

\\Hx-Kx\l\\H 2 -K 2 \\ <e(S) 

where e(5) = E(l/5)5 1 ^ 6 and the function E(x) grows more slowly than any power of x. 

Proof. The only difference here is the requirement that K\ and K 2 be contractions, and the 
construction in (18] does, in fact, produce contractions. 

It is important that the function E(x) does not depend on the dimension n of the matrices. 

□ 

The disk is an easy to understand closed subset of the square, so we expect an easy 
conversion of the quantitative almost commuting Hermitian contractions result to a quan- 
titative result about almost normal contractions. This is foreshadowed by Osborne's result 
about a "bent square" in [36J. (Osborne's result applies to unitaries that correspond to 
^-representations of a subset of the torus homeomorphic to a square.) 

Theorem 2.2. Suppose X is a matrix that is a 5 -representation of the disk. There exists 
X' that is an exact representation of the disk with with 

P"-*'|| < e(5) 

where e(5) = E(1/S)S 1 ^ 6 and the function E(x) grows more slowly than any power of x. 

Proof. This is a corollary of Theorem 12.11 and any improvement on the bounds there will lead 
to an improvement of the bounds here. Let e(S) be the function denoted e(S) in Theorem l2.2l 
Given a contraction X with ||[X*,X]|| < 5, we consider its Hermitian and ant i- Hermitian 
parts. These commute to within 5/2 so there are commuting Hermitian contractions H' and 
K' with 



(X + X*) - H' 



-X + X*) - K' 



< 1(6/2). 



We set X = H' + iK' and X' = f [X I for 



/(*) 



z when \z\ < 1 

A when \z\ > 1. 

\z\ 1 1 
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It is clear that X' is a normal contraction. We easily estimate 

C ||X|| + X -X < 1 + 21(5/2). 
By the spectral mapping theorem 

X'-X < 2e(5/2) 

and so 

||X'-X|| < 4e(5/2). 



□ 



Theorem 2.3. Suppose X is a matrix that is a 5 -representation of the annulus. There exists 
X' that is an exact representation of the annulus with 

||X-X'|| < e{6) 

where e(8) = E(l/5)5 1 / 6 and the function E(x) grows more slowly than any power of x. 

Proof. Let e{$) be the function denoted e(S) in Theorem 12.21 Given X with ||X|| < 1 and 
||X -1 || < 2 and ||[X*,X]|| < 5, we know there is X that is a normal contraction with 



Recall ([231 page 177]) 



So long as 



X -X 
X-X 



X-X 



X" 



< m 



X" 



< 2 X" 



X-X 





we have 








x- 1 


< X" 1 + 


x- 1 - x- 1 


< 2 + 8 


X-X 













This puts the spectrum of X inside an annulus with inner radius (2 + 8e{5)y L . We define 
X' = /(X) for appropriate / so that | < |X'| < 1 and 



|X'-X||<-e(5). 



□ 



We can convert from the annulus to the cylinder rather easily. The spaces are same, but 
the defining relations are different. 

Theorem 2.4. Suppose (U, K) is a 5 -representation of the cylinder by matrices. There exists 
a pair of matrices (U', K') that is an exact representation of the cylinder with 

\\U -U'\\, \\K-K'\\ < e(8) 

where t(5) = E{l/5)5 1 ^ and the function E(x) grows more slowly than any power of x. 
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Proof. Let e(5) be the function denoted e(S) in Theorem 12.31 

Given a unitary U and a Hermitian contraction K, with \\[U, K]\\ < 5, we can form 



Clearly ± < |X| < 1 and 



16 16 



<^\m,u}\\ 



We apply Theorem 12.31 and obtain normal Y with | < |Y| < 1 and 

||F-X|| < e(55/16). 

We convert back with the polar decomposition. Set U' = Y \ Y\~ l and K' — 4 \ Y \ — 31. The 
bounds on |Y| immediately give us — 1 < K' < 1. Since [/' is a unitary that commutes with 
\Y\ it commutes with X'. As to the perturbation estimates, we see 



\K'-K\ 



By [28] we have 



\U'-U\\ < 



= 4|||Y| - |X||| 

= 4||t/'*r-[/*x|| 

< -C/|| +e(5<5/16). 



so 



and 



\U' - C/|| < 3e(55/16) 



\K' -K\\ < 41(55/16). 



□ 



3. Matrices that Almost Represent the Sphere 

Recall that almost commuting triples of Hermitian matrices where the sum of squares is 
almost /, we are regarding as an approximate representation of the sphere. We are still 
studying a two dimensional problem, but have now a non-trivial cohomology class to make 
life interesting. 

3.1. The Index, and when it Vanishes. Useful notation here are the unitless Pauli spin 
matrices, 



"0 


f 




' i 




"1 " 


1 





, 0- 2 = 


-i 


, ^3 = 


-1 



and some indicator functions defined on the real line, 




1 if x > 7, 
otherwise. 



s 
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Definition 3.1. For any triple (Hi, H 2 , H 3 ) of Hermitians we define first another Hermitian 

B (Hi, H 2 , H 3 ) = h + ^J2 a r®H r 



i + h 3 Hi- m 2 

Hi +iH 2 I- H A 



If the H r are n-by-n matrices and | is not in the spectrum of B (Hi, H 2 , H 3 ) then the 
Bott index of this triple is the number of eigenvalues (counted according to multiplicity) of 
B (Hi, H 2 , H 3 ) that are greater than |, minus n, so 

I 0" 




bottOHi, H 2 , H 3 ) = Tr (/i (B (H x , H 2 , F 3 ))) - Tr 



The Hermitian B (Hi, H 2 , H 3 ) is almost idempotent. The Kq groups for C*-algebras are 
generally defined in terms of projections, while for rings one uses idempotents. For a general 
theory of approximate representation of surfaces, the preferred description for indices is in 
terms of approximate projections. 

For this special case of the sphere, another formula demands attention. Let 



S(Hi,H 2 ,H 3 ) = ^a r ®H r 



so that 



H 3 Hi - iH 2 
Hi + iH 2 -H 3 



bott^, H 2 , H 3 ) = ^Tr (/ (S (Hi, H 2 , H 3 ))) . 



We next see that is the triple is a ^-representation of the sphere then 5 < j is enough to 
ensure the Bott index is defined, and as S gets smaller the gap at \ in the spectrum of the 
approximate projection grows larger roughly proportionally. For S(H\, H 2 , H 3 ) the gap is at 



zero. 



Lemma 3.2. If (Hi, H 2 , H 3 ) is a 5 -representation of the sphere for 5 < \ and S = S(Hi, H 2 , H 3 ) 
then US' 2 - III < 45 and 



a(S)C -VI +48, -Vl -45 U VI - 45, Vl + 45 
Proof. From 

S 2 = I <gi (Hi + Hi + Hi) + a 3 ® i [Hi, H 2 \ +a x ®i [H 2 , H 3 ] + a 2 ®i [H 3 , Hi 
we obtain the estimate 

||S 2 -/|| <45. 

The spectral mapping theorem tells us 



o-(S) c{xe 



\x 



1| < 45} 



□ 



The Bott index has appeared in many forms, under different names, in many papers such 
as [7J [13l [HI [311 E21 [33J, [M] . See [9] for a survey of related results in operator theory. 

The Bott index is clearly invariant under conjugation by a unitary. We will see it is very 
stable, is additive with respect to direct sums, and it vanishes when the triple commutes. 
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Example 3.3. For the matrices Hi,H 2 ,H 3 of example (11.21) . bott(ifi, H 2 , H 3 ) 
now show. The approximate projector B is equal to 

1 



1, as we 



(3.1) 
which is equal to 
(3.2) 



B(H U H 2 ,H 3 ) = -I 



Or 



B 



1 



1 



1 



2 2^/S{S+l) 

This is recognizable as the Hamiltonian describing an S"?7(2)-invariant coupling between this 
spin S and an additional spin-1/2. Since B commutes with total spin, there are 2(S + 1/2) + 
1 = 2S + 2 eigenvectors with spin S + 1/2 and 2(S — 1/2) + 1 = 2S eigenvectors with spin 
S — 1/2. The respective eigenvalues are 

{S ± l/2)(5 + 1 ± 1/2) - S(S + 1) - 3/4 



(3.3) 



1/2 + 
1/2 + 



±(S + 1/2) -1/2 



2y/S{S+l) 
« 1/2 + 1/2. 

The difference in the number of eigenvectors with given spin is the index in this case. 

Lemma 3.4. If (Hi, H 2 , H 3 ) is a 5 -representation of the sphere byn-by-n matrices and5 < \ 
then 

bott(# 1; H 2 , H 3 ) - ^Tr (Hi [H 2 , H 3 ]) < 32n5 2 . 

Therefore, if n5 2 < 1/64, bott(#i, H 2 , H 3 ) = Rnd(|Tr [H x [H 2 , H 3 ])), where Rnd(...) 
means round to the nearest integer. 

Proof. Let p(x) = (1/2) (— x 3 + 3x) . If Q is a matrix with spectrum within 7 of ±1 (for 
7 < 1) then p(Q) will have spectrum within 27 2 of ±1. Lemma [3.21 tells us 

(3.4) a (S {Hr, H 2 , H 3 )) C [-1 + 26, -1 - 45] U [1 - A5, 1 + 26} 

Considering the maximum possible errors on the 2n eigenvalues we conclude 

||Tr (/o (S (H u H 2 , H 3 ))) - Tr (p (5 (tf 1; tf 2 , tf 3 ))) || < (2n) (2(45) 2 ) . 

or 



(3.5) 



bott^, H 2 , H 3 ) - ^Tr (p (5 # 2 , # 3 ))) 



< 32n<5 2 . 



Clearly the trace of S(H\,H 2 ,H 3 ) is zero. For the trace of the third power we can drop 
all terms in the product that have two or three indices equal since the trace of any of the 
Pauli spin matrices is zero, so 

Tr (S(H 1} H 2 , H 3 ) 3 ) = a <°s°t ® H r H s H t 

r,s,t distinct 

Tr (±il) Tr (H r H s H t ) 

r,s,t distinct 
= GiTr(H 1 [H 2 ,H 3 \). 
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□ 

Lemma 3.5. Suppose (Hi, H2, H 3 ) and (Ki, K 2 , K 3 ) are triples of Hermitian n-by-n matri- 
ces. Suppose (Hi, H 2 , H 3 ) is a 5 -representation of the sphere with S < \. If 



\\Hi - Ki|| + \\H 2 - K 2 \\ + \\H 3 - K 3 \\ < 

then the Bott index of (Ki, K 2 , K 3 ) is defined and 

bott(Ki, K 2 , K 3 ) = bott(if 1; H 2 , H 3 ). 

Proof. Let S(0) = S (Hi, H 2 , H 3 ) and S(l) = S (Hi, H 2 , H 3 ) . Let 

7 = \\H X - K x \\ + \\H 2 -K 2 \\ + \\H 3 -K 3 \\. 

Clearly \\S — S'\\ < 7. Consider the continuous path S(t) = tS(l) + (1 — t)S(0) and notice 
\\S(t) — 5(0) || < 7. So long as 7 < Vl — AS the gap at zero in the spectrum at zero must 
persist for all t and the indices must be equal. □ 

Lemma 3.6. Suppose 8 < \ and (Hi, H 2 , H 3 ) and (Ki, K 2 , K 3 ) are 5 -representations of the 
sphere by matrices. Then 



bott 



#3 
K 3 



Hi 1 I" H 2 
K x \ ' [ K 2 

= bott(#!, H 2 , H 3 ) + bott(K 1 , K 2 , K 3 ). 

Lemma 3.7. Replacing any one of the H r by —H r flips the sign of the index. 

Lemma 3.8. Suppose Hi, H2, H 3 are three Hermitian matrices so that the Bott index is 
defined. If the H r pairwise commute then 

bott(Hi, H 2 ,H 3 ) = 0. 

Proof. The index is invariant under conjugation by a unitary, we may assume the H r are 
diagonal. The index is additive for direct sums, we may assume n = 1. For real scalars a, b 
and c the matrix 

1 1 + a b + ic 

2 b — ic 1 — a 

has eigenvalues 



-±-Va 2 + b 2 + c 2 . 

We have one eigenvalue above | so the index in this simple case is zero. □ 

This index can vanish for another reason. Lemmas ( 13.9113. 121) consider two cases when the 
index vanishes. Later in section (jGj) we discusses physical motivation for considering these 
cases. 

Lemma 3.9. Suppose S < | and (Hi, H 2 , H 3 ) is a 5 -representation of the sphere. If the Hj 
are real matrices, then 

bott(#i, H 2 , H 3 ) = 0. 
This is a special case of the following. 
Lemma 3.10. Suppose 5 < | and (Hi, H2, H 3 ) is a 5 -representation of the sphere. Then 

bott(H^,H^,Hl) = -hott(Hi, H 2 ,H 3 ). 
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Proof. Consider the unitary 



U 



/ 
-/ 



Then 



Hi 



Hi 



iHl 



Hi 



U * + 2 



I + H! 

H 2 - iH 3 



H 2 + iH 3 
I -Hi 



1 T 


'10' 




/ 



The number of eigenvalues near 1 for the two summands must sum to 2n. Since unitary 
equivalence and transpose preserve all the necessarily-real eigenvalues, we are done. □ 



Definition 3.11. A matrix A is said to be self-dual if ZA T Z = —A, where Z is the block 
matrix 

' o r 
-/ o. 



(3.6) 



Lemma 3.12. Suppose 5 < \ and (Hi, H 2 , H 3 ) is a 5 -representation of the sphere. If the 
H r are self-dual matrices then 

bott(Hi, H 2 ,H 3 ) = 0. 

Proof. This is a special case of the following. □ 

Lemma 3.13. Suppose 5 < \ and (Hi, H 2 , H 3 ) is a 5 -representation of the sphere by 2N- 
by-2N matrices. Then 

bott(#!, if 2 , H 3 ) = -hott(-ZHfZ, -ZHlZ, -ZHjZ). 

Proof. As proven above, bott (ff i, H 2 , H 3 ) = — bott(ff^, ffj, ff J). So, it suffices to prove 
that bo\X(-ZHiZ, -ZH 2 Z, -ZH 3 Z) = bott(ffi, ff 2 , ff 3 ). However, 



(3.7) 



I - ZH t Z -ZH 2 Z - iZH 3 Z 
-ZH 2 Z + iZH 3 Z I + ZH X Z 



Z 
Z 



H 2 - %H 3 
-I + ffi 



= U* 



-I -Hi 

H 2 + %H 3 

I + Hi H 2 + iH 3 
H 2 -iH 3 I- Hi 

I + Hi H 2 + iH 3 
H 2 -iH 3 I- Hi 



Z 
Z 



Z 
Z 

Z 

z 



u, 



where U is the unitary matrix 
(3.8) U = 

Since unitary equivalence preserves the real eigenvalues, we are done. 



Z 

z 



□ 



3.2. Cylindrical to Spherical Coordinates. A key result, in the next subsection, is that 
when the index vanishes for matrices almost representing the sphere they are near matrices 
that exactly represent the sphere. The proof involves a "change of coordinates" into spherical 
coordinates. 

In this subsection we consider the easier change from cylindrical to spherical. 
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Lemma 3.14. Suppose 5 > 0. If (U, K) is a 5 -representation of the cylinder by n-by-n 
matrices then (Hi, H 2 , H 3 ) is a S -representation of the sphere for 

Hi = K, 

#2 = 2 (U^T^K + VI - K 2 U*^ , 
H 3 = % - (-Uy/I - K 2 + VI - K 2 U*^j . 
Proof. These matrices are evidently self-adjoint, and 



Therefore 
and so 
Therefore 



H 2 + iH 3 = UVI - K 2 . 
(H 2 + tH 3 )*(H 2 + tH 3 ) = I-K 2 
Hi + (H 2 + iH 3 )* (H 2 + iH 3 ) = I. 

\\Hl + Hl + Hl-l\\ = \\[H 2 ,H 3 \\\. 



As to the commutators, 

\\[Hi,H 2 + iH 3 ] 



KUVl - K 2 - UVI - K 2 K 
< \\KU -UK\\ 



and 



Since 
we see 
Since 
we have 



(H 2 + iH 3 )* , (H 2 + iH 3 )\ || = ||(1- K 2 ) -U(l- K 2 ) U* \ 

= \\K 2 U - UKK 2 \\ 
< 2 \\KU -UK\\ . 

[(H 2 + iH 3 f , (H 2 + iH 3 )} = 2i [H 2 , H 3 ] 

\\[H 2 ,H 3 ]\\<\\[U,K]\\. 
[Hi, H 2 + iH 3 \ = [Hi, H 2 \ +t[Hi, H 3 \ 



H^^^l + i^,^]!! < \\[u,k\\\, 

and considering real and imaginary parts, we have the weaker estimates 

(3-9) ||[// 1 ,// 2 ]||<||[C/,K]||, 

and 

(3-10) IH^^HII^^II. 



□ 
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3.3. Spherical to Cylindrical Coordinates. 

Lemma 3.15. Suppose (H\,H2,H^) is a 5 -representation of the sphere by matrices. If 
bott(i?i, H2, H3) = then there is a unitary U so that (U, H3) is a \/~85 + 25 -representation 
of the cylinder and 



(3.11) 



U (I — Hf) 5 — (ifx + iH 2 ) 



< 2^2r/ + 2r/. 



Proof. Let P = B(Hi, H2, H3). The index vanishing means there is a unitary W so that 



W* 



D x 
D 



for Di a diagonal matrix within 25 of I. Therefore 

/ 



P-W* 







W 



w 



< 25. 



' A 


B ' 




" / 




















" A 


B ' 




" A 


B ' 




" / 


" 

























Define A and B by 



Then 
(3.12) 
and 

so 

(3.13) 

From (13.121) we get one exact relation, 

AA* + BB* 

and from (13.131) several approximate relations, 



W. 



' A 


B ' 


* 


" A 


B ' 


= W 


' I 


" 























w* 





'AO' 


* 


'AO' 






B 




B 


-P 



< 25. 



1/ + Ij7 3 
2 2 3 



2 2 



i/ - iff, 

2 2 



< 25, 

< 25, 

< 25. 



(3.14) A* A 

(3.15) A*B - 

(3.16) B*B 
In particular, 

\\A*A + B*B - I\\ <45. 

Recall we can insist on a unitary in the polar decomposition of a matrix, although it may 

not be unique. Also notice that X = U (X*X)^ implies X = (XX*) * U. See §83 in [17J, for 
example. Thus there are unitaries Z and V so that 

A = Z(A*AY = (AA*Y 2 Z 



14 

and 

Next we find 
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B = V{B*B)i = (BB*)*V. 



V[h + ^H 3 j V* - AA* 

\ ' ( --I + ^H 3 + B'D ) \ ' - \ 'B BY - , L V . I 



< 



-I - -H 3 - B*B 
2 2 

-/ - -H 3 - B*B 

2 2 6 



\\VB*BV* + AA* - 
\\BB* + AA* - I\\ 



so 



(3.17) 
and 

















V 






- aa* 




Z*AA*Z 


-( 






a*a- 

















so 



(3.18) 

These two equations now tell us 



Z A.-VZ - ( -I + -H 3 
2 2 



< 25. 



\\(Z*V)*H 3 (Z*V)-H 3 \\ 

= \\Z(I + H 3 )Z*-V(I + H 3 )v*\\ 

< \\Z (I + Hi) Z* - 2AA*\\ + \\2AA* - v (I + 

< 85. 

Let [7 be the unitary U = Z*V. We just showed (U,H 3 ) is an 85-representation of the 
cylinder. From equation (I3.17P we get 

V {^1 + ^H 3 y V* - {AA*f <V25, 

c.f. [37] or [3]. To equation (13.161) we apply equation (1) in [3] to produce the estimate 



(B*BY 



h - 1 -h, 

2 2 s 
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and now 



U{l-HlY - (H x + iH 2 



U (k I+1 2 H > 



2 2 6 



< 



-H + -H 2 
2 2 



u(~I + ~hX (B*Bf 



Z*V (-1+ -H 3 ) 2 (B*B)^ - Z* (AA*Y V (B*B)^ 



A*B - ( -Hi + -H 2 

1 2 -r 2 



< 



+ 



h - iIh, 

2 2 s 



2 2 3 



(B*B) 
- (AA*)^V 



A*B 
<V85 + 25. 



-Hi + -H 2 
2 2 



□ 



Theorem 3.16. Suppose (Hi, H 2 , H 3 ) is a 5 -representation of the sphere by matrices. If 

bott(Hi, H 2 , H 3 ) = 0, 
then there are commuting Hermitian matrices H'^H'^H'^ with 

(H[) 2 + (H' 2 f + (Htf = I 

and 

\\K-H r \\ <e(5) 

for all r, where e(5) = E(l/5)5 1 ^ 12 and the function E(x) grows more slowly than any power 
ofx. 

Proof. Let e(S) be the function denoted e(S) from Theorem 12. 41 
By Lemma [3.151 there is a unitary U so that 



\u*mu-H,\\ < 65 



and 



U(I- Hi) 5 - (Hi + iH 2 ) < 2^65 + 65. 
Theorem 12.41 produces unitary V and Hermitian contraction K that commute and with 

\\V-U\\ < e(65) 

and 

\\K-H 3 \\ <e(65). 



16 
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Define H[, H 2 , H' 3 by 

H[ = - (VVl - K 2 + Vl - K 2 V*} = Re (yy/l - K 2 ^ , 
H' 2 = % - (-Wl-K 2 + Vl-K 2 V*^j = Im (Vy/l - K 2 ^j ., 



H' 3 = K. 



Lemma [3. 141 implies that V and K commute and satisfy 



Of course 

Finally 
\\Hi + iH 2 - H[ - iH' 2 \ 



m 2 + (H' 2 f + (h>: 



|#1 -#i|| <e(65). 



Hi + iH 2 - VVl - K 2 



< 



H 1 + iH 2 -UJl-H 2 



UJl-H 2 -Wl-K< 



< 2^65 + 65 + \\U - V\ 



I -HI- Vl-K 2 



< 2V65 + 65 + e (65) H 
< 2^65 + 65 + e (65) + ^~2\\H 3 -K 



HI) ~ (1 



K 2 )\\ 



< 2V65 + 65 + I (65) + y/2e (65). 



□ 



Remark 3.17. The power of 1/12 in Theorem 13.161 can be improved upon. It is possible to 
modify the construction in [18] to improve the approximation of one of the Hermitians as 
the cost of weakening the approximation of the other. This leads to an asymmetric version 
of Theorem 12.41 and an improvement to Theorem 13.161 

Theorem 3.18. On Hilbert space HI there are bounded Hermitian operators Hi, H 2 and H 3 
and finite rank projections Pi < P 2 < . . . so that: 

(1) the strong limit of the P n is the identity I; 

(2) the H r commute; 

(3) forr = 1,2,3 we have lim \\[P n , H r }\\ = 0; 

ra— >oo 

(4) if K n i, K n>2 and K n3 are commuting Hermitian operators then 

\\K nA - P n HiP n \\ + \\K n . 2 - P n H 2 P n \\ + \\K n ,z - P n H 3 P n \\ -> 1. 



Proof. The idea is to put the matrices H n ^, H n2 , H n ^, from Example 11.21 down the diagonal 
of an infinite matrix, but "doubling" as in [22]. Our almost commuting projections will cut 
one of the double blocks in half. 
The index of 



H n ,i 






-H Ui i 



H n ,2 







-H n2 



H n ,3 







H n3 
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is zero, so we can approximate these by X n , Y n and Z n that are commuting 2n-by-2n Her- 
mitians whose squares sum to one and with 

H nil 
-H nA 



X n 



0. 



Y 



H n , 2 

-# n ,2 

H n , 3 

-H n , 3 



0. 
0. 



Let 



Q, 



I 




Let X, Y, Z and P correspond to the block diagonal matrices formed out of 

Xi, X2, . . . , X n -x, X n , X n+ i, X n+ 2, . . . 
Yi, Yii • • • , Y n _i, Y n , Y n+ i, Y n+ 2-, . . . 
Z\, Z-ii ■ ■ ■ , Z n _i, Z n: Z n+ i, Z n+ 2, ■ ■ ■ 



and 

We have 



h, ■ ■ ■ , hn-2, Qni 0, 0, ... . 



\XP n — P n X\ 



X„ 



I 




X r , 







and similarly ||yP n — ^^11 ~ > an d \\ZP n — P n Z\ 



I 


0. However, the index of 



P n X P n) PnY Pni PnZP n 



is 1, so cannot be approximated by commuting Hermitians. 

4. Matrices that Almost Represent the Torus 



□ 



Definition 4.1. Suppose U and V are unitaries and ||Z7V — VU\\ < 2. The winding number 
invariant u(U, V) of (U, V) is the winding number of the closed path in C \ {0} given by the 
formula 

t ^ det ((1 - t)UV - tVU) . 

Let log(-) denote the branch of the logarithm continuous except on the negative reals. Let 
Tr denote the trace on M n , normalized by Tr(J) = n. The following is due to Exel ( [T21 p. 
213]). 

Lemma 4.2. So long as \\[U,V]\\ < 2, 



u{U,V) 



2tt 



Tr (-i\og(VUV*U*)) 



Proof. There is a homotopy between two paths from UV to VU. The first path is the linear 
path. The second path sends t to e ^g(vuv*u*) UV . See [EH p. 213] for details. □ 

The following should be obvious. 

Lemma 4.3. If U and V are commuting n-by-n unitaries then co(U, V) = 0. 

This invariant is very stable. 
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Lemma 4.4. IfUo,Ui and Vq,Vi are n-by-n unitaries with 

max(||[[/ ,V ]|| J[^i,K]||)<<5 



then 
implies 



|t/o-t/i| 



|Vb — Vill < 2-5 



u(U ,V )=u(U 1 ,V 1 ). 
Vo — Vi\\ . Consider the following, defined for (s,t) in the unit 



Proof. Let r\ — ||£/o — Ui\\ 
square, 

X s>t = s t U Vo - s tiV Uo + SitoUiVi - sitiVit/i, 

where t = 1 — t, t± = t, sq = 1 — s and Si = s. We want this to be invertible. The "corners" 
are unitaries, so if this stays within distance 1 of a corner it will be invertible. Standard 
estimates show 

\\X s , t - XJ^W < s Q tt \\[U , V }\\ + 8l tt \\[U U K]|| + 8j+1 (\\Ui ~ U \\ + ||K - KID 

\\X,, t - VjUiW < s t \\[U , V ]\\ + Sl t HfC/i, Vj]\\ + s j+1 (Hf/x - U \\ + \\Vl - V \\) 
where the j + 1 is to be performed mod-2. Putting in the assumptions on the norms we find 

\\X Sjt - UjVj\\ < hS + s j+1 r) 

\\X S)t - VjUjW < t S + s j+1 r] 
so there is always one corner to which the distance is less than ^5 + \r\. 



□ 



Example 4.5. Let Uj = e " 3 and 



L0 2 



S n 





1 





1 



Then Q, n S n Q,^S* equals e » I so u>(S n ,Q n ) = 1. 

39] , but we have a large lower bound on the distance 



This is an old example 
to commuting unitaries. 

Theorem 4.6. If U and V are commuting n-by-n unitaries then 

"K 



\U-Sr, 



\V-n n \\ > 2 1 - sin 



n 



Proof. 



\u-s n \\ + \\v-n n \\ <2-\\[s n ,n r 



2iri 

e n 



2 - 2 sin ( - ) . 

.n, 



□ 
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This index shares a lot of properties with the Bott index from j 
with respect to conjugation by a unitary. Also 

u(V, U) = oj(-U, V) = u(U, -V) = -oj(U, V) 

and 

n ' " 1 w(U u V x )+u,(U 2 ,V 2 ). 



It is obviously invariant 













( 




> 




) 



In fact, this is equal to an invariant directly based on i^-theory. 

There are no "really nice" maps from the torus to the sphere, but there are smooth maps 
that are one-to-one over most points of the sphere. For present purposes, the smoothness is 
not so important. 

On map from the standard torus in C 2 to the unit sphere in M 3 that is one-to-one over 
most points of the sphere, and is piecewise smooth, is 



,2mdi 



) = (/i(0i),<7i (6i) + h (£i)cos 



h (00 sin (0 2 )) 



where 



if x < \ 
ifs<±' 



9xK x ) 




ifs<± 



hi(x) 







ifx<h 



if x < \ 



V-8 + 24x - 16x 2 
1, and from here it is easy to check that T does take 



Notice gxh\ = and /f + g[ + h{ 
values on the unit sphere. 

The restriction of V to the open set determined by (1/2) < 9\ < 1 gives a bijection onto 
the sphere minus the poles. Points in the complement are sent by T to a half-equator. 

This function give us a means to manufacture an approximate representation of the sphere 
out of an approximate representation of the torus. Unfortunately, the norm of a commutator 
|| [X, f(Y)]\\ depends rather poorly on the norm of \\[X, Y]\\ when X and Y are normal and 
all we know about / is that is is piece-wise linear or smooth. 

Suppose U and V are approximately commuting. If we write out the approximate projec- 
tive B based on the three almost commuting Hermitians (in some order), it looks like 

f(V) g(V)+h(V)U~ 
g(V)+U*h(V) I-f(V) 

where /, g and h are the continuous, real-valued functions on the circle defined by 



Q(U,V) 



g ( e 2m9j _ ig^ ^ ^ an( j ^ (e 2m6 ^j — lfi x (#) . The three Hermitians that almost represent the 

/ + 2f(V), 2Re (g(V) + iH{V)) , 2Irn (g(V) + iH{V)) . 

Lemma 4.7. There is a 8\ > so that for all unitary matrices U and U with \\[U, V]\\ < 5±, 
the Hermitian matrix Q(U, V) does not have | it its spectrum. 



-,2-KiO 

sphere are 
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Definition 4.8. For all unitary matrices U and U with \\[U, V]\\ < S\, define k{U, V) as —n 

plus the number of eigenvalues, counted with multiplicity, of Q(U, V) that are greater than 
1 

2' 

Notice that Q(U, V) varies continuously in U and V. There is a variation on the index that 
lacks this feature, but it has cleaner formulas. 

In [TJ] it was determined that there is a So, again unspecified, so that || [U, V] \\ < S implies 

yfe(v) - (t(V)) 2 U 
U*Jl(V) - (l(V)f £(V) 

has spectrum that does not contain |, where £(e 2mx ) = x for x G [0, 1). We define K\(U, V) 
as —n plus the number of eigenvalues, counted with multiplicity, of Qi(U, V) that are greater 
than The following is proven in [14]. Taking into account Lemma l4~2l we have four methods 
to compute the index. 

Theorem 4.9. There is a 63 > so that for all unitary matrices U and V with || [U, V]\\ < S3, 
we have equality of the indices, 

w(U,V) = k(U,V) = K!(U,V). 

We have no quantitative version of the following, and no proof of the following that does 
not utilize the theory of C*-algebras. We refer readers seeking a proof to [151 Corollary M3] 
and [TUl Theorem 6.15]. 

Theorem 4.10. For every positive e, there is a positive S less than 2 so that \\[U,V}\\ < S 
and uj{U, V) = for unitary matrices implies that there exists commuting unitary matrices 
Ui and V\ with \\U — U\\\ < e and \\V — Vi\\ < e. 

5. Applications 

5.1. Lattice Problems and Definition of Wannier Functions. We describe an appli- 
cation of the previous results to a physics problem: whether or not there exist localized 
Wannier functions for a two-dimensional insulator. We begin by describing the physical sys- 
tem, and the appropriate mathematical description of this problem. We then connect the 
existence of Wannier functions to the absence of an obstruction to approximating almost 
commuting matrices. See also [26J. 

We consider a system of non-interacting fermions moving in a tight binding model with the 
sites arranged on the surface of a sphere (this topology is chosen to give a two-dimensional 
system with no boundary) and with short-range hopping terms of bounded strength. 

Since the fermions are non-interacting, the Hamiltonian of the system is 

(5.1) H = J2*lH lj y 3 , 

where are fermion creation and annihilation operators on site i. The matrix if is a 

Hermitian matrix. Since we consider insulating systems, we assume that the Hamiltonian 
has a spectral gap AE, so that, without loss of generality, all eigenvalues are either less than 
or equal to —AE/2 or are greater than or equal to AE/2. We let P denote the projector 
onto the space spanned by the eigenvectors with eigenvalues less than or equal to —AE/2. 



Qi (u, v) = 
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A basis of Wannier functions means an orthogonal set of states which spans the range of 
P. An example of a basis of Wannier functions would simply be the set of eigenvectors with 
eigenvalues less than or equal to —AE/2. A local basis of Wannier functions means a basis 
of Wannier functions in which some locality requirements are imposed on the basis vectors: 
for example, given a metric on the lattice, for each vector most of the norm of the vector 
should be concentrated on sites near some given site. Usually, the eigenvectors of H will 
not be local in this sense. Below, we give a few precise, but slightly different, mathematical 
definitions of "localized Wannier function". 

To describe the fact that the interactions are short-range, we need to describe where each 
site is located on the surface of the sphere. To do this, we introduce matrices X, Y, Z which 
are diagonal matrices. For each site i, the corresponding diagonal matrix elements Xu, Yu, Zu 
denote the x, y, z positions of that site on the surface of a sphere. We impose the condition 

(5.2) X 2 + Y 2 + Z 2 = L 2 , 

where L is the radius of the sphere. We define the distance between any two sites by: 



(5.3) dist(i,j) = J(Xii - X n f + (Y u - Y 3J f + (Zu - Z 



\2 



03 1 



Define a distance between a site i and a triple of coordinates x = (x, y, z) by 



(5.4) dist(z,£)) ee ^ ( Xll - x) 2 + {% - y) 2 + {Z u - zf. 

To mathematically describe the assumption of short-range interactions, we assume that 
Hij = for dist(z, j) > R for some range R. To describe the assumption of bounded strength 
interactions, we assume that \\H\\ < J for some interaction strength J. The conditions are 
sufficient to imply a Lieb- Robinson bound for the dynamics [29] [T9| [35] . These conditions 
imply that 

(5.5) \\[X,H}\\ < 2v LR , 

for vlr = RJ, and similar bounds for ||[Y,if]||, \\[Z, H]\\. The subscript LR refers to Lieb- 
Robinson; this velocity that we define here can be shown to be an upper bound on the ve- 
locity of propagation of excitations in according with the usual definition of a Lieb-Robinson 
velocity. 

We will be interested in the case where L » vlr/AE below. We now derive a bound on 
|| [P, X] || , || [P, Y] || , || [P, Z] || . Because of the spectral gap, we can write 

(5.6) P = AE J dtf(AEt)exp(iHt), 

for any function f(t) such that the Fourier transform, f(u>) obeys f(ui) = 1 for u < — 1/2 
and f(u>) = for u > 1/2. Then, 

(5.7) \\[X,P)\\ = AE\\[X,Jdtf(AEt)exp(iHt))\\ 

< AE J dt\f{AEt)\\\[X,exp{iHt))\\ 

< 2AE J dt\f{AEt)\v LR \t\. 
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We now choose any specific such f{ui) which is at least twice times differentiable for all u, so 
that f(t) decays faster than 1/t 3 for large t. Then, the integral on the last line of Eq. (15.71) 
converges and 

(5.8) || [X, P] || < const, x v LR /AE, 

for some numeric constant, and the same bound holds for \\[Y, P]\\ and ||[Z,P]||. 
We now define matrices H r by 

(5.9) Hi = PXP/L, 

H 2 = PYP/L, 
H 3 = PZP/L. 

We now show that these matrices approximately represent the sphere: 

Lemma 5.1. The matrices Hi, H 2 , H3 defined in Eq. ( 15. 9\) form a 5 -representation of the 
sphere with 

(5.10) 5 = const, x (v LR /LAE) 2 . 

Proof. Define X n = (1 - P)X(1 - P),X 12 = (1 - P)XP,X 2 i = PX(l - P), and X 22 = 
PXP = Hi. Define Y\\ = (1 — P)Y(1 — P), and so on, so that we can write 

(5.11) X = (i 11 * 12 ) , Y f Yn Yl2 



\X 2 i X 22 J \Y 2 i Y 22/ 

Then, P[X,Y]P = X 21 Y 12 - Y 21 X 12 + [X 22 ,Y 22 \. However, since [X,Y] = 0, this means 
that [X 22 ,Y 22 ] = -X 21 Y 12 + Y 21 X 12 , so ||[#i,# 2 ]|| < ||^2i|| \\Y 12 \\ + ||r 2 i|| ||X 12 ||. Note that 
||^2i|| < ||[-X",-P]|| < const, x v LR /AE, and similarly for ||X 12 ||, ||F 2 i|| || Vi 2 ||- So, \\[H U H 2 }\\ < 
const, x (v LR /LAE) 2 . 

Similar bounds hold for the commutators || [H 2 , H 3 ] || , || [H 3 , Hi] \\ . 

Finally, 

\\H 2 + H 2 + H 2 -I\\ = \\PX(1-P)XP + PY(1-P)YP + PZ(1-P)ZP\\ 
(5.12) < \\px(i -p)xp\\ + ||py(i -P)YP\\ + \\PZ(1 -P)ZP\\ 

< ||(1 - P)XP\\ 2 + 11(1 - P)FP|| 2 + ||(1 - P)ZP|| 2 

< ||[P,X]|| 2 +||[P,f]|| 2 + ||[P,Z]|| 2 

< const, x (v LR /LAE) 2 . 

□ 

From now on, we work in the subspace spanned by eigenvectors with eigenvalues less than 
or equal to —AE/2; i.e., the subspace onto which P projects. In this subspace, H r are a set 
of almost commuting Hermitian operators that almost square to unity. The question now 
is: do there exist a set of localized Wannier functions? We give three possible definitions 
of this, and relate these definitions to the ability to approximate H r by exactly commuting 
matrices. 

Definition 5.2. A set of exponentially localized Wannier functions with localization 
length £ is a set of orthonormal vectors, {v a }, spanning the subspace onto which P projects, 
such that for each vector v a the following property holds. Let v® denote the coefficient of v a 
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in basis element i. Let x a = (v a ,Xv a ), y a = (v a ,Yv a ), and z a = (v a ,Zv a ). We require that, 
for all a and all D, 

(5.13) Yl \vt\<eM-D/0- 

i,dist(i,(x a ,y a ,z a )>D 

Definition 5.3. A set of Wannier functions localized to length l[ oc is a set of orthonor- 
mal vectors, {v a }, spanning the subspace onto which P projects, such that the following prop- 
erty holds. Let x a = (v a ,Xv a ), y a = (v a ,Yv a ), and z a = (v a ,Zv a ). Let x a = (x a ,y a ,z a )). 

Let X' = £ a z |i; )(t; |, let Y' = £ a z |u )(u°|, and let & = Ea^KX^I; We require that 
for any w in the subspace onto which P projects that 

(5.14) \(X-X'))w\ <hoc\w\, 

\(Y-Y'))w\ <l loc \w\, 
\(Z - Z'))w\ < h oc \w\, 



Definition 5.4. A set of Wannier functions weakly localized to length l\ oc is a set of 
orthonormal vectors, {v a }, spanning the subspace onto which P projects, such that for each 
vector v a the following property holds: 

(5. 15) (v a , X V) - {v a , Xv a ) 2 < lf oc , 

(v a ,Y 2 v a )-(v a ,Yv a ) 2 <ll c , 
(v a , Z 2 v a ) - K, Zv a f < l 2 oc . 

Definition (15.31) implies definition (15.41) ; to see this note that Eq. (15.141) for w = v a im- 
plies Eq. (15.151) . Under one assumption, definition (15. 2p implies definition (15.31) . This is an 
assumption about the number of points in the original lattice, as in the following lemma; 
this assumption expresses the two-dimensionality of the original problem. This next lemma 
unfortunately is fairly tedious in the details, given the simplicity of the resulting estimate. 

Lemma 5.5. Suppose that, for any x = (x,y,z) and any I, the number of sites j with 
(x — Xa) 2 + [y — Ya) 2 + [z — Za) 2 < I 2 is bounded by C\ + c 2 l 2 , for some constants c 1; c 2 . 
Then, 

(5.16) \(X - X'))w\ 2 < const, x £ 2 [ Cl + c 2 £ 2 ] 2 
and similarly for Y —Y' and Z — Z' . 

Proof. The assumption on the number of sites in the lattice implies a bound on the number 
of vectors v a with dist(a;, x a ) < I as follows. Define 

(5.17) P(x,l)= K)^!- 

a,dist(x,a; a )<Z 



\vf\ 2 
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Then, for any D, 

(5.18) Tr(P(f,Z)) 

dist(x,x a )<l 

= E E k*i 2 

i a 

dist(z,i)<Z+D dist(z,£ a )<Z dist(:r,i)>Z+Z) dist(x,x a )</ 

= E E ki 2 + E E 

i a 2 a 

dist(x,i)>l-\-D dist(x,x a )>D 

< c 1+ c 2 (i+Df+ J2 E ki 2 

i a 

< ci + c 2 (Z + P) 2 + Tr(P(f, Z)) exp(-D/0- 

Picking P = £, we see that Tr(P(af, Z), which is the number of such vectors v a with 
dist(x, x a ) < I, is bounded by const, x (ci + c 2 (Z + £) 2 ). Next, for w = J2 a A(a)v a , we 
have 

(5.19) \(X-X')w\ 2 

= ~M°)A(b) (v a , (X - x a ) (X - x b )v h ) 

a^b 

= Ma)A{b) KvK*ii ~ * a ) (X* - * b ) 

aj^b i 

< \ A (a)A(b)\ MrfWiXti - x a )(X u - x b )\ 

aj^b i 

< |A(a)A(6)| ^exp(-dist(z,x a )/e) exp(-dist(^)/£)l(^ - x a ){X u - x h )\. 

aj^b i 

The sum over i in the last line of the above Equation (I5.19P is equal to 

(5.20) exp(-dist(«,f a )/Oexp(-dist(«,f 6 )/OI(X ii -x a )(X ii -x f, )| 

i,dist(i,x a )<2dist(x a ,x b ) 

+ exp(-dist( i ,f a )/Oexp(-dist(«,f 6 )/OI(X ii -x a )(X ii -a; fe )| 

i,dist(i,x a )>2dist(i a ,x b ) 

< const, x (exp(-dist(f a ,f b )/Odist(f a ,f 6 ) 2 [ci + c 2 (dist(f a ,f fc )+0 2 ] 

ofc ^ dist(i,£ a ) ^^fc+l 
OO dist(x a ,^^) — 



+ E E exp(-dist(z,f a )/e)exp(-dist( i ,f fc )/e)|(X, i -x a )(X J ,-x 6 )|) 

fe=l i 

< const, x (exp(-dist(f a ,f b )/e)dist(f a ,f 6 ) 2 [ci + c 2 (dist(f a ,f fc ) +0 2 ] 

OO 

x(l + ^exp(-2 fc dist(x°,x b )/£)(2 • 2 fe ) 4 )) 

k=l 

< const, x exp(-dist(x a , f 6 )/0(£ + dist(x a , f 6 )) 2 [ci + ^(dist^, f 6 ) + f) 2 ] 
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Combining (15.191) with (I5.20p gives 
(5.21) 

\{X-X')w\ 2 < const. x^|A(a)A(6)|exp(-dist(x^f b )/0(^+dist(x a ,f b )) 2 [ci+c 2 (dist(f a ,f 6 )+0 2 ] 

a^b 

Combining this with (15.181) . we have: 
\(X-X'))w\ 2 

< const, x ^|A(a)A(6)|exp(-dist(f a ,f 6 )/0(^ + dist(f a ,f 6 )) 2 [ci + C2(dist(f a ,f b )+0 2 ] 

a^b 

< const, x ^[(|A(a)| 2 + |A(6)| 2 )/2] exp(-dist(f a , x 6 )/^)^ + dist(f a , f 6 )) 2 ^ + ^(dist^, f 6 ) + 2 ] 
= const, x ^|A(a)| 2 ^exp(-dist(x a ,x b )/0(^ + dist(x a ,f b )) 2 [c 1 + c 2 (dist(x a ,f b )+0 2 ] 

a byta 

< const, x £ 2 [ci + c 2 £ 2 ] 2 

□ 

5.2. Obstructions to Wannier Functions and Almost Commuting Matrices. We 

claim that definition (I5.3P is equivalent to the ability to approximate Hi,H 2 , H 3 by exactly 
commuting matrices H[,H! 2 ,H! i . Consider first the direction of the implication that (15.31) 
implies the ability to approximate by exactly commuting matrices: simply set 

(5.22) H[ = ^x a \v a )(v a l 

a 

H' 2 = Y,y a \v a )(v% 

a 

H' 3 = ^z a \v a )(v a \. 

a 

Then, \\Hi~ H[\\ = m&x W) \ w \ =1 \(Hi—H[)w\ = m&x w ^ w \ = i\(PXP— X')w\/L = max u , ) | lu | = i|(P.X'— 
X')w\/L < max^, i w i=i \(X - X')w\/L < l/L. Similar bounds follow for \\H 2 - H' 2 \\ and 
11-^3 — LL 3 \\. So, Eq. (15.141) implies the ability to approximate by exactly commuting matrices 
up to error l/L. To see the converse implication, that the ability to approximate by exactly 
commuting matrices implies definition (15.31) . let the vectors v a be basis vectors in a basis in 
which H' r are exactly diagonal. Then, for any w, 

(5.23) \(X-X')w\ = \]\{PXP - X'))w\ 2 + |(1 - P)Xw\ 2 

= L^\{H l -H[)w\ 2 + \[P,X]w\ 2 
< - H[\\ 2 +4(v LR /AE) 2 \w\. 

Combining these implications, the presence of an index obstruction to approximating by 
exactly commuting matrices implies an obstruction to definition (15.31) . which implies (under 
the assumption above about the number of points in the original lattice) an obstruction to 
finding exponentially localized Wannier functions. Conversely, the ability to approximate 
H r by exactly commuting matrices implies the ability to find Wannier functions obeying 
definitions fl5.3U5.4j) . 
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Combining lemma (15. ip with theorem (13.161) . the absence of an index obstruction im- 
plies the ability to find Wannier functions localized to length Le(const x (vlr/ LAE) 2 ) = 
L(vlr/ LAE) 1 / 6 £7 (const, x LAE /vlr)- This length scale is asymptotically smaller than L. 
We leave as an open problem the question of whether the absence of an index obstruction 
implies the ability to find Wannier functions localized to length vlr/AE. We also leave as 
an open problem the question of whether the absence of an index obstruction implies the 
ability to find exponentially localized Wannier functions; such a result is known in the trans- 
lationally invariant case[5]. Note in this regard that under the assumptions of finite range 
and interaction strength and spectral gap AE, it is possible to prove that the matrix ele- 
ments of P are exponentially decaying: \Pij\ < exp(— dist(z, where £' is proportional 
to Vlr/AE. This proof uses standard techniques to prove locality of correlation functions in 
gapped systems [20l [21] and is based on using a smoother function f(t) in Eq. f!5.6p . 

The index obstruction can be computed for several examples. In an ordinary band in- 
sulator, the obstruction vanishes. In numerical results below, we give applications to a 
lattice realization of a quantum Hall system on the surface of the sphere where the index is 
non- vanishing. 

5.3. Numerical Simulations. We have performed numerical simulations to illustrate the 
usefulness of this index obstruction for studying physical systems such as a quantum Hall 
effect on the sphere. The index we consider has the advantage, compared to more usual Chern 
number obstructions calculated on a momentum torus[5], that it does not require translation 
invariance, and it also does not require averaging over parameters of the Hamiltonian on a 
flux toruspQ, as in the Chern number calculation of the Hall conductance. Compared to the 
noncommutative geometry approach discussed in j2] , we have a method that will adapt to a 
variety of surfaces and is applicable to finite size systems. 

We have considered the following system. The choices that we made are deliberately 
somewhat arbitrary: we wanted a modest sized system describing free particles moving 
on the surface of a sphere in the presence of a roughly uniform magnetic field exiting the 
sphere, but we wanted to illustrate the robustness of this index even in a system with no 
carefully chosen symmetry. We considered a total of 560 sites on the surface of a sphere. 
The sites were distributed on 29 different latitudes, such that all sites in a given latitude 
had the same angle from the north pole (and hence had the same z coordinate). The 
angles 9 describing the latitudes were evenly spaced from 7r/30, 27r/30, ...,297r/30. On each 
latitude, the number of different sites equal to the floor of a constant times cos(#), for 
some constant giving 560 total sites, with the angles <fi of the sites evenly spaced from 
to 27T. The matrices X, Y, Z were chosen to equal the x, y, z coordinates of each site, with 
Xu = sm(0j) sin(0j), Y i{ = cos(6 l i ) sin(<^), Z i{ = cos(0;). 

The Hamiltonian H was chosen so that Hij = if the distance between sites i and 

j, measured as J (Xu — Xjj) 2 + (Yu — Yjj) 2 + (Zu — Zjj) 2 , was greater than a maximum 

range, which we chose to be V0.07 ~ 0.26. Otherwise, the matrix element was chosen 
to equal —Jijexp(iujij), where Jy was the strength of the interaction and uoij was a phase. 
We set Jij equal to —1 plus a constant 5 times a random number chosen independently for 
each pair i, j and uniformly between —0.5 and 0.5. The phase was chosen to mimic the 
effect of a magnetic field. We picked 

(5.24) Uij = n monopo i e * ((pi - <f)j) * cos((6»i + 0j)/2), 
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Figure 5.1. Spectrum of the Hamiltonian for 5 = (black line) and 5 = 1 
(red line). The right column is a detail of the figure near the band gap. 

where n monopo i e is an integer describing the net flux leaving the sphere. We chose n monopo i e = 
100 to make the net flux slightly smaller than the number of sites. 

Diagonalization of the Hamiltonian for 5 = 0, 1 gave the energy spectra shown in Fig. (15. 3p . 
Note that for 5 = there are very few eigenvalues between roughly —2.3 and —0.3 and there 
are no eigenvalues between —1.3994 and —0.7296. It is important to understand that the 
term "band gap" can be used in two different ways in physics. One way means that there 
is an interval (E min , E max ) containing no eigenvalues. This is referred to as a "strict band 
gap." The other use is that there is an interval containing very few eigenvalues, with the 
corresponding eigenvectors being localized. Such eigenvectors are referred to as "mid-gap 
states." If there is a strict gap in the spectrum of H, of order unity, and if the range of H is 
much less than unity (for example, 0.07) in our case, then the gap can be used to prove that 
P approximately commutes with X, Y, Z as discussed above. In the example we considered, 
however, there is not a very large strict band gap: if we want to consider that the energy 1 
lies in the middle of a strict band gap, then all we can say is that there are are no eigenvalues 
in the interval (—1.3994,-0.7296). However, even without a large strict band gap, if the 
mid-gap states are indeed localized, then the projector P will still approximately commute 
with X, Y, Z and so the commutators of PXP, PYP, PZP will still be small. We will see 
numerically below that this is the case for our problem. 

There are 200 eigenvalues less than —1, so the projector P onto states with energy less than 
— 1 has rank 200. We computed the spectrum of B for matrices PXP, PYP, PZP. The index 
was equal to unity, so that there were 201 eigenvalues close to unity and 199 close to zero. 
The smallest eigenvalue close to unity was equal to 0.9763 and the largest eigenvalue close 
to zero was equal to 0.0146, so that there is a very clear separation between the eigenvalues 
close to zero and those close to unity. The largest commutator was [PYP, PZP] , with 
|| [PYP, PZP] || « 0.0298021. Thus, the matrices are very close to commuting as claimed. In 
Fig. (15. 3p we plot the eigenvalues of B. 
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Eigenvalue number Eigenvalue number 

Figure 5.2. Spectrum of B for 5 = 0, projecting onto energies above 1. 
The left and right columns plot the eigenvalues close to unity and to zero, 
respectively. There are 201 eigenvalues on the left and 199 on the right. 

The case with 6 = 1 was similar. There were, in this case, 201 eigenvalues less than — 1 
(so one eigenvalue crossed unity at some value of 5 between and 1). However, the index 
remained equal to unity, and the smallest eigenvalue of the B close to unity was 0.9642 
while the largest eigenvalue close to zero was 0.0242. The largest commutator was again 
[PYP, PZP], with || [PYP, PZP] || « 0.036728. Thus, the matrices are very still close to 
commuting. 

We can explicitly check that the index is insensitive to small changes in the energy as long 
as we do not enter the band of delocalized states. Projecting instead onto eigenvalues less 
than —2, there are a total of 195 states, the index is still equal to unity, and the eigenvalues 
of B which were closest to 0.5 were 0.305 and 0.9562. Projecting onto eigenvalues less than 
—3, there are a total of 163 states, the index is still equal to unity, and the eigenvalues of B 
which were closest to 0.5 were 0.3571 and 0.6123. Projecting onto eigenvalues less than —4, 
there are a total of 90 states, the index is now equal to zero, and the eigenvalues of B which 
were closest to 0.5 were 0.4555 and 0.6657. Thus, as expected, when we move away from 
the band gap, the index ceases to be well-defined since the matrices cease to approximately 
commute. 

Finally, we can check that we approximately represent the sphere, namely that PXP 2 + 
PYP 2 + PZP 2 is close to the identity in the subspace projected onto by P. The smallest 
eigenvalue was 0.934833 when projecting onto states with energy less than —1 and the 85-th 
largest eigenvalue was still greater than 0.99, but the smallest eigenvalue was only 0.447221 
when projecting onto states with energy less than —4. 

5.4. Relation to Hall Conductance. The index calculated numerically for the lattice Hall 
system above is clearly closely related to the Hall conductance. The formula in lemma (13.41) 
expresses an approximation to the index in terms of a trace Tr (PXP [PYP, PZP]). Consider 
a family of Hamiltonians H with increasing L, with a uniform lower bound on the spectral 
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gap, and uniform upper bound on v LR . Suppose that the dimension of the Hamiltonians, 
n, is proportional to L 2 , as is natural for a two-dimensional system. Then, by lemma (15. ip . 
the matrices H r form a 5- representation of the sphere, with nS 2 oc 1/L 2 . This means that, 
by lemma (13.41) . Tr(PXP[PYP, PZP\) is within 1/L 2 of an integer. That is, this trace is 
approximately quantized. 

This trace is closely related to the Kubo linear response formula for the Hall effect: con- 
sider applying an electric potential to the sphere which varies uniformly between the north 
and south poles. This amounts to adding a term ^ Z^> to the Hamiltonian. If the Hall 
conductance is positive, this will drive a current in the counterclockwise direction around 
the sphere (and in the clockwise direction for negative Hall conductance). Thus, at pos- 
itive x-coordinate, one would expect to see a current in the positive y-direction, and at 
negative x-coordinate one would expect to see a current in the negative ^-direction. If the 
Hamiltonian H is proportional to a projector, then this response is proportional to the trace 
Tr (PXP [PYP, PZP]), up to numeric constants, and factors of the electric charge. Thus, 
we prove approximate quantization of the Hall conductance in spherical geometry for non- 
interacting electrons for such Hamiltonians which are proportional to projectors. Perhaps 
with more work it will be possible in this way to prove quantization of the Hall conductance 
in spherical geometry for non-interacting electrons for arbitrary gapped Hamiltonians. 

6. Matrices With Additional Reality Constraints 

In this section, we consider further the case in which the matrices H r are assumed to be 
either real or self-dual. By the results above, since the index vanishes in this case, it is 
possible to approximate almost commuting real or self-dual matrices by exactly commuting 
matrices. However, we can ask a further question: is it possible to approximate real or 
self-dual matrices by exactly commuting real or self-dual matrices? 

We begin with physical motivation for considering this problem. Based on the physical 
intuition, it is natural to conjecture that there is a Z2 obstruction to approximating almost 
commuting self-dual matrices by exactly commuting self-dual matrices, and that there are 
no other constructions. We then verify one of these conjectures: we construct this Z 2 index, 
prove that if the index is nontrivial then there is a lower bound on the distance to exactly 
commuting self-dual matrices, and construct an example with a nontrivial index. We then 
finish with precise statements of our other conjectures. 

6.1. Physical Motivation. In case the Hamiltonian Ti, has time reversal symmetry, the 
matrix P will have the same symmetry. The possible cases of interest correspond to different 
universality classes in random matrix theories. We discuss three classes here, corresponding 
to the GUE, GOE, and GSE classes. For previous application of these universality classes 
to classifying different insulating phases of free fermions, see [38], in particular table II. 

In the GUE case, H, has no time reversal symmetry, P is a Hermitian projector, and H r 
are Hermitian matrices with no further symmetry constraints. In the GOE case, Ti has time 
reversal symmetry, and the time reversal symmetry operator squares to unity: this describes 
a spin-0 particle in a time reversal symmetric situation, or a spin-1/2 particle with time 
reversal symmetry and no spin-orbit coupling. In this case, H r are real, symmetric matrices. 
In the GSE case, 7i has time reversal symmetry, and the time reversal symmetry operator 
squares to minus one. This describes a spin-1/2 particle with strong spin-orbit coupling. In 
this case, the H r are self-dual. 
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In the GSE case, the index vanishes due to the time reversal symmetry as explained 
above. There is therefore no index obstruction to approximating H r by exactly commuting 
H' r . However, it is natural to look for a Z 2 obstruction to approximating three almost 
commuting H r by self-dual H' rl because we know that there is a Z 2 index characterizing 
different translationally invariant phases of free fermions with symplectic symmetry in two 
dimensions [25], [23] - Such topologically nontrivial phases are expected, by their nature, to be 
stable to perturbations of the Hamiltonian which break translational symmetry, as discussed 
in [38] and have been observed experimentally in HgTe/(Hg,Ce)Te quantum wells[27j. 

6.2. Index Obstruction. We construct the index by finding a unitary transformation that 
makes S(Hi, H 2 , H 3 ) anti-symmetric and then taking the sign of the Pfaffian of this matrix. 
From a physical point of view, the existence of this unitary transformation is not surprising: 
the self-dual operation can be regarded as a time-reversal symmetry operation, and a similar 
time-reversal symmetry can be applied to the a matrices used to construct S(H 1 , H 2 , H 3 ). 
Under these combined time reversal symmetries, S(H\, H 2 , H3) changes sign; however, since 
there are two spin- 1 /2s, the time reversal symmetry operator squares to unity and hence, up 
to a basis change, is equivalent to transposition. We now show this. 

Lemma 6.1. Let H r be self-dual. Define the matrix B(Hi, H 2 , H 3 ) by 



Proof. Note that for any r, aj = —cr 2 a r a 2 , while Hj = —ZH r Z, by the assumption of 
self-duality. Also, U T = U. Thus, 



(6.1) 



B(H U H 2 , H 3 ) = U*S(H U H 2 , H 3 )U, 



where the unitary U is defined by 



(6.2) 




Then, the matrix B is anti-symmetric. 



(6.3) 




r 




r 




r 




r 



Note that (Z g> o 2 f 



—I, so 



(6.4) 




-U*, 
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and 
(6.5) 

Thus, 
(6.6) 



(Z®a 2 )u* = -^={1 + Z <g> a 2 )U 



V2' 
U. 



U(z ®a 2 ^{^H r ®a r ^{z ®a 2 ^jU* = -U*(j2 H r®°r)u 

r r 

= — B(Hi, H 2 , H 3 ). 



□ 



Definition 6.2. We define the index bott(ifi, H 2 , H 3 ) for self-dual matrices H r by 
(6.7) Qt(ifi, H 2 , H 3 ) = sgn(Pf ff 2 , ff 3 ))), 

where Pf is the Pfaffian and sgn(x) = 1 for x > and sgn(:r) = —1 for x < 0. If 
Pf # 2 , # 3 )) = 0, the index bott (if i, if 2 , #3) is not defined. 

Lemma 6.3. Consider any continuous path of self- dual matrices, H r (s), where s is a real 
number, < s < 1. Suppose that for all s, the matrix B(Hi, H 2 , if 3) /10s non-vanishing 
determinant. Then, boTt(ifi(0), if 2 (0), H 3 (0)) = bott(ifi(l), ff 2 (l), ff 3 (l))- 

Proof. The determinant of B(Hi, if 2 , if 3 ) is equal to Pf (B(Hi, if 2 , H 3 )) 2 . As long as the 
determinant does not vanish, the Pfaffian does not vanish and hence cannot change sign. □ 

Lemma 6.4. If Hi, H 2 , if 3 are self-dual and exactly commuting and bott(ifi, if 2 , if 3) is 
defined, then bott (if 1, if 2 , if 3) = 1. 

Proof. Note that bott is multiplicative under direct sum of matrices. Also, bott is invariant 
under symplectic transformation of the if r (these transformations preserve the property of 
being self-dual). Finally, for any commuting self-dual matrices if r , we can find a symplectic 
transformation which makes the if r diagonal; the diagonal entries of these matrices come 
in pairs which are equal. That is, if the if r are 2n dimensional matrices, then after this 
symplectic transformation then each if r is equal to the direct sum of n different 2-by-2 
matrices which are proportional to the identity. Thus, it suffices to consider the case in 
which the H r are real scalar multiples of the 2-by-2 identity matrix. If if r = a r I then 



(6.8) 



B(Hi, H 2 , H 3 ) 



Z ® <T 2 ) ar ® a r) ( J + Z 

ia 3 Z ® 0-3 + ia 2 I ® a 2 + ia 3 Z ® a\ 

( %ol\ ia 2 ia 3 \ 

—ia\ — ia 3 ia 2 

ia 2 ia 3 — ia\ 

\— ia 3 ia 2 %ol\ 



°2 



This matrix has Pfaffian equal to 

(iai)(-iai) - (ia 2 )(ia 2 ) + (ia 3 )(-ia 3 ) 



2 2 

a*! + a 2 



a 3 . 



□ 
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We now show that for any 5 > 0, there exist self-dual matrices H r which form a 8- 
representation of the sphere with bott(ifi, H 2 , H 3 ) = — 1. This example essentially consists 
of two copies of the matrices in (11.21) . with opposite Bott index between the two copies. 

Example 6.5. Let n = 2S + 1. Consider the spin matrices S 1 , S 2 , S 3 for a quantum spin S 
with Hi = g> S 1 /^S(S + 1), H 2 = a ( 2 l) <8> S 2 /,/S(S + l), H 3 = I« ® S 3 /y/3{3+ 1). 
We will be using the a matrices in two different ways in this example: first, to form the 2n 
dimensional matrices H r from the n dimensional matrices S r , second to form the matrix B. 
We use erW to refer to the first case, and to refer to the second. Note that Z = —io~^®I- 
We use 7^ to refer to the 2-by-2 identity matrix. The matrix S 2 is anti-symmetric, while 
S 1 , S 3 are symmetric. 
It is easy to see that 

\\H r , H s \\ < 1/ S, 

so that the matrices form a ^-representation of the sphere with 5 = 1/S. 
We claim that bott(#i, H 2 ,H 3 ) = -1. 

Proof. We now compute the Pfaflian. Since the index depends only on the sign of the Pfaflian, 
we ignore constant factors which are real and positive. We have 

(6.9) B(H U H 2 ,H 3 ) 

= const, x (I - Z ® a ( 2 } )B(H u H 2 , H 3 )(I + Z® a ( 2 ] ) 

= const, x a™ <g> (s x ® a { 3 ] + S y ® a { 2 } - S z <g> cr{ 2) ) . 

Unitarily conjugate this matrix B(Hi, H 2 , H 3 ) by the orthogonal transformation (1/V2) J (g> 
(/( 2 ) +?o"2 2 ' 1 ), giving the matrix const, x cr^ S r (Tr. Since this orthogonal transformation 

has determinant +1, the Pfaflian is unchanged. 

The matrix a 2 ® Er S r °~r 2 ^ is anti-symmetric and equals 



3.10) 



iEr^r 



The Pfaflian of this matrix is equal to the determinant of i^^Vr, which is equal to 
i 4S+2 times the determinant of E r S T a r . The matrix E r S^cv has 2(S + 1 /2) + 1 = 2S + 2 
positive eigenvalues and 2(S — 1/2) + 1 = 2S negative eigenvalues, as computed in example 
( 11.2p . Thus, the sign of the Pfaflian is equal to 

(6.11) ,4S+2 ( _ )2 S = _ L 

□ 

Lemma 6.6. Suppose (Hi, H 2 , H 3 ) and (Ki,K 2 ,K 3 ) are triples of self- dual, Hermitian n- 
by-n matrices and suppose (Hi,H 2 ,H 3 ) is a 5 -representation of the sphere with 5 < 1/4. 
If 

\\Hi - Ki\\ + \\H 2 - K 2 \\ + \\H 3 - K 3 \\ < VT^45 

then 

boTt(K 1; K 2 , K 3 ) = bdft(#i, H 2 , H z ) 
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Proof. We showed in the proof of Lemma [3.51 that the line segment from S(H\, H 2 , H 3 ) to 
S(K\, K 2 , K3) passes through invertibles. It follows that the line segment from B(H\, H 2 , H 3 ) 
to B(Ki, K 2 , K3) passes through skew-symmetric invertibles and so by lemma (16.31) the Pfaf- 
fian does not change sign. □ 

As a corollary, the distance in operator norm from the matrices in Example ( 16. 51) to the 
nearest exactly commuting triple of self-dual matrices is at least yfl — A/S. 

6.3. Conjectures. We believe that the absence of this index obstruction implies that is 
possible to approximate almost commuting self-dual matrices by exactly commuting self- 
dual matrices. Before addressing this issue, we need to understand the two matrix case in 
the presence of additional symmetry. We thus raise the following conjectures which generalize 
Lin's theorem: 

Conjecture 1. For all e > 0, there exists a 5 > such that, given any real, symmetric 
matrices A, B with ||L4, B}\\ < 5 and ||A||, ||S|| < 1, there exist real, symmetric matrices 
A', B f , with [A', B'] = and \\A - A'\\, \\B - B'\\ < e. 

Conjecture 2. For all e > 0, there exists a 5 > such that, given any self-dual, Hermitian 
matrices A, B with \\[A, B]\\ < 5 and \\A\\, \\B\\ < 1, there exist self-dual, Hermitian matrices 
A f ,B', with [A 1 , B'} = and \\A - A'\\,\\B - B'\\ < e. 

The conjectures we make regarding index obstructions are: 

Conjecture 3. For all e > 0, there exists a 5 > such that, given any real, symmetric 
matrices Fi,if 2 ,#3 with ||[#i,# 2 ]||, ||[#2, # 3 ]||, Il[#3,#i]|| < and 

(6.12) H* + Hl + Hi = I, 

there exist real, symmetric commuting matrices H[, H' 2l H' 3 , with ||-f^2 — H' 2 \\, ||-f^3 — 

"3 II < ^ 

Conjecture 4. For all e > 0, there exists a 5 > such that, given any self-dual Hermitian 
matrices H U H 2 ,H 3 with bott(# 1; H 2 , H 3 ) = 1 and \\[H U H 2 ]\\,\\[H 2 , H 3 }\\,\\[H 3 , < 5, 
and 

(6.13) Hl + Hj + Hl = I, 

there exist self-dual Hermitian commuting matrices H[,H 2 ,H^, with 11^ — H[\\, \\H 2 — 
H 2 \\, \\H 3 - H' 3 \\ < e. 

7. Discussion 

We have given quantitative error bounds on the ability to approximate three almost com- 
muting matrices by three exactly commuting matrices under the assumption of a vanishing 
index assumption. We have related the ability to approximate these matrices to the ability 
to find localized Wannier functions in a physical system, and we have demonstrated that it 
is readily possible to numerically calculate this index for such systems. 

We have constructed a Z 2 index obstruction to approximation of almost commuting self- 
dual matrices by exactly commuting self-dual matrices, and raised additional conjectures 
regarding almost commuting matrices in the case of real C*-algebras. Finally, we note that 
table II of [38] lists 10 different classes of matrices and different index possibilities in various 
dimensions. We believe more generally that each of the obstructions in d = 2 in this table will 
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correspond to a particular obstruction to approximating almost commuting matrices; i.e., 
we expect that for almost commuting matrices in class Dili there will be a Z 2 obstruction 
to approximating them by exactly commuting matrices of class DHL Such obstruction is left 
for future work. 
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